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Abstract 

The dynamical spin fluctuations in a two-dimensional square lattice in its para- 
magnetic phase are examined within the framework of Random Phase Approxima- 
tion(RPA). Itinerant carriers with spin interact with each other via an antiferromag- 
netic Heisenberg interaction. 

Then there appear three fundamental scattering processes; a) scattering between 
parallel spins, b) scattering between antiparallel spins and c) scattering with spin-flip. 
To examine how these scattering processes affect the dynamical spin fluctuations, we 
pick up carefully all possible combination of RPA diagrams in a consistent manner and 
take the spin rotational symmetry into account. Then it becomes clear that we have to 
take up a sequence of the irreducible single loop which in itself is modified due to the 
particle-hole ladder type vertex correction. We set up the Bethe-Salpeter equation for 
the vertex correction and show that this can be solved in a closed form due to separable 
nature of the antiferromagnetic interaction. We evaluated numerically the effect of the 
vertex correction and found that the correction is negligibly small. 

Therefore we propose that in an itinerant Heisenberg system, including the t-J 
model as its derivertive, the simplified RPA, where the irreducible single loop is un- 
renormalized, works very well. This conclusion strongly support the simplified treat- 
ment which is widely used in High-T c problem. Moreover the present formalism enables 
us to proceed further microscopic calculations on the magnetic properties in the current 
High-T c problem. 



* kishine@tansei.cc.u-tokyo.ac.jp 
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1 Introduction 



The purpose here is to explore the dynamical spin fluctuations in a paramagnetic state of 
a two-dimensional square lattice electron system in terms of the itinerant Heisenberg model 
within the framework of the naive random phase approximation (RPA). In this model the free 
motion of the itinerant carrier, electron or hole with spin, is modified by scattering processes 
via an antiferromagnetic Heisenberg interaction^. Then there appear three fundamental 
scattering processes; a) scattering between carriers with parallel spins, b) scattering between 
carriers with antiparallel spins and c) scattering between carriers with spin-flip. These 
processes affect each other and appear in the renormalized spin fluctuating processes in all 
possible ways. So far, however, little attention has been given to the naive treatment of this 
problem. 

The similar problem has been studied, in connection with current High-T c problems^], 
in terms of the t-J model which can be regarded as the derivative of the present model. 
In the t-J model where the effect of non-double occupancy condition due to the strong 
correlation is incorporated and consequently the magnetic processes are governed by the 
spinon's degree of freedom. Then there exist scattering processes between the free spinons 
via an antiferromagnetic Heisenberg interaction. In this sence the same problem occurs 
concerning the spin fluctuating processes as in the case of the present model. 

Indeed, within the framework of a mean field picture of the t-J model, the spin fluc- 
tuations have been intensively studied by Tanamoto, Kohno, and Fukuyama (T.K.F)[|J 
where the dynamical spin-spin correlation function, Xt.k.f(<z) = Xspmon[l + ^qXspmont^)] -1 , 
plays the central role. Here Xspinon(g) is the irreducible single loop of a free spinon and 
J q = 2J(cosq x a + cos q y a) is an antiferromagnetic Heisenberg interaction in a momentum 
space. 

However the form of the correlation function Xt.k.f{q) include only one single scattering 
channel and, from diagrammatic point of view, consists of a sequence of the irreducible single 
loop which is unrenormalized. Moreover this form of the susceptibility doesn't contain the 
spin indexies as x al3 ■ The simple form similar to Xt.k.f((z) has been widely used in the 
context of the High-T c problem. 

The present work is motivated by a naive question; "when we consider all the possible 
scattering processes within RPA, what kind of the scattering processes contribute to the 
dynamical spin fluctuations and consequently the simplified structure of Xt.k.f(<?) should be 
modified or not?" . 

In the present paper we put aside the effect of the strong correlation and concentrate 
only on the structure of possible scattering processes which produce the dynamical spin 
fluctuations within the framework of RPA. 

1 Although the present work was motivated by the usual treatment of the t-J model. We thought it better 
to aboid calling the present model "t-J" model, since the present model neglect the non-double occupancy 
condition due to the strong crrelation. Our standpoint is similar to RefJlJ in which the spinon loop is 
regarded as ordinary hole and the effect of the non-double occupancy condition is not taken into account. 
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Generally the naive RPA corresponds to picking up the so-called ring diagram which is 
a sequence of the irreducible single loop. Then the irreducible single loop naturally includes 
the particle-hole exchange scattering processes as the lowest order vertex correction ||. This 
treatment fully reproduces an equation of motion method. 

However, in specific problems, how to choose the proper diagrams is model-dependent. 
For example, in the paramagnon theory based on the Hubbard model@, the transverse spin 
fluctuations consist of one renormalized single loop incorporating the particle-hole exchange 
scattering and the longitudinal spin fluctuations consist of a sequence of the unrenormalized 
single loop. This situation comes from the short range nature of the Hubbard interaction. 

As an another example, in the Coulomb gas problem, the paramagnetic susceptibility 
contains a single loop modified due to the exchange scattering^]. 

It will be shown that in the present case both of the transverse and longitudianl spin 
fluctuations consist of a sequence of the irreducible single loop modified due to the particle- 
hole exchange scattering processes. 

Since we consider a paramagnetic phase, it is also important to pay our attention to the 
spin rotational symmetry of the theory. 

In §2 we present the model. In §3 we review the general formalism for the spin fluctuation 
and construct the naive RPA. Then we show that Bethe-Salpeter equation for the vertex 
correction can be solved in a closed form due to the separable form of the exchange interac- 
tion. We obtain the expression for the spin fluctuations in a closed form, which is different 
from the form where the irreducible single loop is unrenormalized. In §4 we discuss the spin 
rotational symmetry in the present treatment. Finally in §5 we present the numerical results 
and discuss quantitative difference between the present treatment and the simlified one. 

As a result we propose that the simplified treatment for the spin fluctuations of a square 
lattice which has been widely used works fairly well in the present framework. 



2 Model 

We start from the itinerant Heisenberg model, 

H = Ho + Hint, (1) 

Here 

= E&cLck-r. (2) 

kcr 

is the kinetic Hamiltonian where c-^J (ck CT ) is a creation (annihilation) operator of an itinerant 
carrier with the momentum k and the spin projection a. 
Further 

H% n t = 2 J^t/qSq ■ S_ q , (3) 



3 



represents the antiferromagnetic Heisenberg interaction between the nearest neighbor spins 
where 



J q = J (cos q x a + cosg y a), (4) 

denotes the momentum dependent spin-spin interaction. The spin fluctuation operator with 
momentum q is defined by 

S q = 9 C l.+c l ,a a 'oil3Ckl3- (5) 
k,a,3 

where a denotes the usual Pauli matrixes and we set H = 1. 

Bearing the high-T c problem in mind, we include not only the nearest neighbor hopping, 
t, but also the next nearest neighbor hopping, t', and then we have 

6c = £k — A* = — 2t(cos k x a + cos k y a — a cos k x a cos k y a) — /x, (6) 

where a is the lattice constant and a = —2t'/t. Here < a < 1 is a parameter which 
characterlizes the geometry of the Fermi contour. Furthermore /i denotes the chemical 
potential. 

The dispersion (H) produces the density of states (DOS); 



V a (e) = ^ ._ K 



{a-e/2t)\ 



7T 2 v^T^i L ^ 2(2 + J ' 

where K(x) denotes the elliptic integral of the first kind. The derivation of (0) will be 
presented in Appendix A. The carrier density in the electron picture is determined by 

n = = / deV a (e)f(e), (8) 

JM J-2t(2-a) 

where f(e) = (e^r 1 + l) -1 . N c and represent the number of carrier and lattice sites 
respectively. Then n = 1 corresponds to the half-filling. 

This type of DOS shows the logarithmic van-Hove singularity at e = —2a for < a < 1 
and the square root singularity for a = 1, which originates from the saddle points located 
at the four points (±7r, 0), (0, ±7r) in the first Brillouin zone. The case of a = which 
corresponds to LSCO compounds leads to the perfect nesting. As the parameter a becomes 
nearer to a = 1, the curvature of the convex along the T-Y line near the saddle point becomes 
flatter. Then a = 1 corresponds to the perfectly flat curvature of the convex. This type of 
the saddle point is called "the extended saddle point" which produces the power singurality 
in DOS pi . In the case of YBCO, a = 0.9 corresponds to the real data|§, which is very near 
the case of the extended saddle point. In Fig. 1. we show the energy contour, corresponding 
DOS, and the structure of the saddle point near the Y point (0, n) for typical values of a. 
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3 Transverse Spin Fluctuation 
3.1 General Formalism 

We consider the magnetic response in a square lattice electron system within the framework 
of the present model. Whole information on spin dynamics is contained in the dynamical 
wave number dependent magnetic susceptibility, 

X ^{q) = |^re^<T T [^(q,r)^(0)]>=iE<, r ^ V (?)< (9) 

where 

r^ ;V (g) = fdre^}_ £ < T T [4, /1 (r) Ck , +q ,,(r)c| c+qiP (0) Ck , A (0)] >, (10) 
J0 iV k,k' 



denotes the spin dependent polarization function where a,/3 = +,—,z and fi,u,X,p =f, | 
and iu n = 2niTn is a bosonic Matsubara frequency. Throughout the present paper q denotes 
a four vector q = (q, iu) n ). 
Furthermore 

cljr) = exp(rW)c k)Q (0) exp(-rW), (11) 

represents imaginary time dependent creation operator where we set H — 1 and Jcb — 1. 
< • • • >= Tr(e~ ■ • .)/Tr^=r- denotes the thermal average under the full Hamiltonian. 



T 

By noting 



^-(S;K-(!!!)--(i-0- (12) 

the longitudinal and transverse spin susceptibility can be obtained as 

X + ~(q) = T^^(q), (13) 

x -+(g) = r ^(?), (14) 

x «(,) = l [r TT;TT (g) + r U;U (g) _ T^(q) -T^(q)] = ^(q) - xHq)\- (15) 

where x^' (<?) = T a ^ a '^'(q). 

Throughout this paper we consider the paramagnetic pahse and therefore the dynamic 
spin susceptibilities must satisfy the rotational symmetry relation in spin space; 

?X ZZ = X + -- (16) 
The dynamical susceptibility of a free carrier system can be given by 

XirUv) = -^Y,Go(k + q)Go(k). (17) 

k 
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Here 

Go(k) = . 1 t , (18) 

l£n - £k 

is a thermal green's function of a free carrier where e n = (2n + l)7rT is a fermionic Matsubara 
frequency. Furthermore from now on the sumation over the four momentum means 

By taking the Matsubara summation we obtain 

Xfree(<?) = E o TTt t \o 5 (20) 

where fermionic Matsubara frequency and 

A k , q (T) = tanh(|±^)-tanh(^), (21) 
denotes the thermal extinction factor. 

3.2 RPA 

Now we consider the dynamical susceptibiliy which is modified due to J q . We calculate the 
dynamical susceptibility within the framework of RPA. 

The fundamental processes induced by J q are produced through 

Sq ■ S_ q = ^ (^q *^-q + ^q ^-q) + ^q^-q- (22) 

By rewriting this expression in a second quantized form, we get three fundamental vertexes 
for the scattering processes between propagating carriers with spin. 



The first term of fl22[) produces the scattering process with spin flip (type-a). The second 
term of (|22f) produces two fundamental processes; scattering processes between parallel spins 
(type-b) and anti-parallel spins (type-c). These processes are shown in Fig. 2. 

Here we first consider the transverse spin fluctuation. The problem of the rotational 
symmetry in spin space will be taken up in the next section. The string type summation as 
shown in Fig. 3. produces the form, 

X + ~(1) = *; ( ?' M . (23) 

where Xo(l) is the irreducible single fermion loop which cannot be cut into two pieces by 
removing an interaction line and which can be written in the form 

UQ) = ~ E 9{k)T(k, k + q)Q{k + q), (24) 

k 
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where T(k, k + q) is the triangle vertex inserted in a single loop. 

Furthermore Q{k) is a dressed green's function determined through the Dyson equation 



w i = t?o(*r i -E(*). ( 25 ) 

where XX^) denotes the self energy part. 

We can see that in Fig. 3. only the type-a processes can contribute to the ring diagram 
for the transverse susceptibility within RPA. Here we should note that, even within the 
framework of RPA, the type-c processes can produce the exchange scattering between a 
electron and hole inside a single loop. Consequently we have to include the particle- hole 
ladder process in the vertex T(k + q,k). These type of diagrams are characteristic to the t-J 
model while in Hubbard model the string type summation cannot appear in the transverse 
fluctuation. 

Therefore we have to include the string type summation and the ladder type summation 
simultaneously to get the naive result within RPA. 

To treat the compatibility of the self-energy and the vertex correction, it is sufficient 
to consider only the Fock term to ensure the Ward-Takahashi identity [[LOU . This procedure 
corresponds to the simplest case of Baym-Kadanoff's conserving approximation^^. Then 
J2{k) can simply be written by 

E(*) = -|^EAk-k')W), (26) 

k' 

which can be reduced to the form 

3 

y^(fe) = —-JX(T) (cos k x a + cos k y a), (27) 



where 



, „ cos k'a + cos k'a 

ie n < - e v + fi - £(#) 



I r dk! 

-J ^p(cosA;> + cosA;;a)tanh[(£ k , - // + E(fc')/2T]. (21 



Here we've performed Matsubara summmation and took into account the T>^ symmetry of 
the system. We can obtain the self energy by solving (p8|) in a self-consistent manner. X(T) 
depends weakly on the temperature T and turns out to have positive value. It follows from 
this situation that the Fock term, as usual, enhances the effective mass. Now the energy 
dispersion of a single carrier is slightly modified from flf^) to the form 

3 

£k = —2[t JX(T)] (cos k x a + cos k y a) + 2a cos k x a cos k y a. (29) 
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The triangle vertex T(k, k + q) satisfies the Bethe-Salpeter equation for a particle-hole 
channel as shown in Fig. 4., 

T(k + q,k) = l + ~J2 J^'S(k')T(k', k' + q)G{k' + q). (30) 

k' 

We can solve (|30|) in the closed form as 

r(fc, k + q) = 1 + JC (q) (cos k x a + cos k y a) + JS(q) (sin /c x a + sin fc^a) (31) 

where 

r( , -|xi(g)[i + \JxM\ + ^[|x 2 (g)][|x 4 (g)] , . 

W " " [i + ^x 3 (g)][i + ^x 5 (g)]-[^X4(g)] 2 ' 1 j 

9r x -\x2{q)[l + jjXaCg)] + A\xM][\xA{q)] , . 

W [i + ^x 3 (g)][i + ^x 5 (g)]-[|Jx4(g)] 2 ' 1 } 

We may leave the details of the derivation of this formula to Appendix B. Here Xi(q)> ' ' *> 
Xb(q) are defined by 



xi(q) = - 




X2(q) = - 


T 


--y 


xs(q) = - 


T 


X4(g) 


T 





y i H h ^ 

^(cos k x a + cos k y a) 

k 

(sin fc^a + sin k y a)Q{k)Q{k + g), 
T 

Xs(<?) = --77 ^(sinA^a + sin k y a) 2 Q{k)Q{k + g). 

Therefore we can obtain the result 

Xo(q) = Xo(q) + JC{q)xM + JS(q)xz(q)- (35) 

where 

xo(q) = ~T,S{k + q)G(k). (36) 

We note here that for the commensurate spin-fluctuation with q = Q = (n /a,7r/a), 
X2 = Xi — due to the _D 4 symmetry of the square lattice, and then consequently S(q) = 0. 

Directly from the above expressions (34), we can expect the correction terms in the 
(|35|) are negligiblly small. Therefore we expect in our scheme the vertex correction coming 
from the exchange scattering processes can be neglected. We will confirm this situation 
numerically in §5. 



S 



4 Rotational Symmetry in Spin Space 



When we consider the spin fluctuations in a paramagnetic phase, to ensure the consistency 
of the theory, the dynamical spin susceptibilities have to satisfy the symmetry relation (|16|) 
in spin space. Here we prove the rotational symmetry of the theory in the present scheme. 
Our concern is to show the relation (O), 



2X 2 



X" 



(37) 



First we show that the spin rotational symmetry is satisfied in a single loop level. Namely 



2x3" 



Xo = Xo, 



(38) 



where xo was given by (p4j) . Then we have to include various intermediate spin configurations 
between a given initital and final spin configuration. We consider a single loop of the n-th 
order term, Xo(n)i with respect to J q . There can appear the type-a and the type-b vertexes 
in all the possible manner, as is shown in Fig. 5. We have to take summation with respect to 
all the possible intermediate spin configurations and pick up a lot of diagrams. Fortunately 
we can obtain the relation between Xo^ (n), Xo( n ) an d Xo(n) by order- by-order consideration. 
The results follows. 
When n = 2m, 



When n = 2m + 1, 



~U 

^0(2m) 



fc=0 
m—1 



~+- 9 2fc+l (~i 

^0(2m) A 2m^2fc+l- 



(39) 
(40) 



fe=0 



~Tt 

Xo(2m+l) 



~1l 

Xo(2m+l) 



"Xo(2m+l) ^ 2k 2m+lC 2 k, 



k=0 
m—1 

~Xo(2m+l) X/ 2 
fc=0 



2fc+l 



2jn+lC2fc+l, 



(41) 
(42) 



where ^C/ = k\/[(k — /)!/!] is a binomial coefficient. Here, for example, we derive the formula 
(|39|). In this case n = 2m and thererfore the single loop can include the even number of the 
type-a processses. Then all the other vertexes are the type-b vertexes. We consider the case 
when there are 2k type-a processses and 2m — 2k type-b vertexes. Since the type-a vertex 
gives the facor J q and the type-b vertex gives the facor J q /2, if we replace all the type-a) 
vertexes simply by the type-b vertexes, there appear the factor 2 2k . On the other hand there 
are 2mp2k ways of locating the type-a vertexes inside a loop. As a result the corresponding 
expression has a factor 2 2k 2 m C 2 k and therefore we obtain d39|). Here XoU) ^ s a sm gl e 1°°P 
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corresponding to the transverse fluctuation of the same order with respect to J q . Now by 
noting that 

m m—l 

Er,2k ri \ ~* o2fc+l n f1 o\2m _ -i 

* 2m ( ~'2k — 2-^i 1 2m<-^2fc+l — \\ ~ ^) — ±, 

fc=0 fc=0 

and 

m m—l 

Eo2fc o2fc+l n — ( ~\ r^m+l _ i 

£ 2m+l ( -'2k — / j ^ 2m+l < ~ / 2fc+l ~ \ l ^ ^) ~ ~ L i 

k=0 k=0 

we can see that for arbitraly order, 

Xo(n) ~ Xo(n) = Xo(n)7 (43) 

and consequently 

Xo - Xo l = Xo~ = Xo- (44) 
Therefore we obtain the result (|3g|) which would be expected. 

Next we consider the string type series for the longitudinal spin fluctuations. In the 
string processes, in this case, there can appear the type-b and type-c vertexes. Then the 
m-th order term with respect to J q can be expressed by 

c™' = (-^) m J2 (-ly^i-iy*-* 4 ■■■(-i) a2m - i - a2k (45) 

Xo ai Xo 2a3 ---Xo 2ma '- (46) 

This process is shown in Fig. 6. 



By using ( p"5|) we can obtain the longitudinal spin susceptibility as 



1 oo 

X ZZ = ^Y.^-C^). (47) 

Z m=0 



First we take the summation over o\ and in ([5]) to get 

^ ff3,"4,-,U2ra 

(xr - xriixr - xo aa3 )xr 5 ■ ■ ■ x?^' 

_ ^_'l3_y n ^_]^°"3-°"4 . . . ^_^C2ro-l- Cam 

^ 0"3,0-4,---,(T2 m 

X0U0 _ Xo JXo Xo 

Here we used the relation (f44"|). 



By repeating the summation procedures over the pair (<7 2 j_i, o"2i), we get the result for 
the z-th procedure as 

^00' _ ^ ^q-jm2«^i+l ^ -Q0 2 i+l-CT 2 i + 2 . . . ^ ]^02m-l-02m 

^ 0"2i+l,0"2i+2,---,0-2m 

(C-Xo^^-xr' (49) 
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Finally we obtain the relation 

C^' ^(-JXXoixf ~Xo a ' a ')- (50) 

Therefore we get 

^ m=o ^ fc=o ^ 1 + J q^o 

which is just one half of the result for the corresponding transverse susceptibility ( P^) . Thus 
we have proven the expected relation ([57]). 



Numerical Results and Conclusion 



In this section we show numerically that the correction terms in (pq) can be neglected 
and justify the anticipation that a simplified treatment where the irreducible single loop is 
unrenormalized works well. 

For this purpose we concentrate on evaluation of the terms JC(q)xi(q) and JS(q)x2(q) 
in (1351). 



From now on, to proceed numerical work, we set the parameters J = t/4. 



l.Self Energy 

First we present the numerical results of the self-consistent equation for the self energy; fl2"8[) . 
We performed the integration over the first Brillouin zone divided into 64 x 64 uniform mesh. 
We set the carrier density n = 0.85. In Fig. 7(a) we show the temperature dependence of 
the factor X(T) in fl28|) for a = and a = 1. We can see that the self energy depends very 
weakly on the temperature over the wide range < T < 200 [K]Q. 

Next in Fig. 7(b) we show the carrier density dependence of X(T). It follows from 
this that in case of a = the effective mass becomes monotonically heavier as the system 
approaches the perfect nesting. On the other hand in case of a = 1 the effective mass has 
the highest value near the carrie density n ~ 0.7 

In any event the Fock term gives rise to the negative correction to the hopping integral t 
of the order |JX(T) ~ 0.005t which can safely be neglected. It is true that we should take 
the Fock term into account to ensure the consistency of the theory, but it becomes negligiblly 
small in magnitude and can be neglected. 



2. Irreducible Single Loop Xo(q) 

Next we are concerned with the correction terms in the single loop. If the single loop 
can well be approximated by the unrenormalized one, Xo(<?)> we can replace the irreducible 

2 Here we put t - 130 [meV] 
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polarization by the polarization for the free carrier system (fL7|), since as shown above the 
self-energy correction can be neglected^. In such a case our treatment reproduces the effective 
RPA in Ref.[|| . 

Here we fix the temperature to T = 0.04t ~ 60[K] and take the corresponding lowest 
matsubara frequency iO\ ~ 32[meV]. The integration is performed over the first Brillouin 
zone divided into 128 x 128 uniform mesh. 

First in Fig. 8 we present the numerical results for the modified polarizations (|34]) in the 
first Brillouin zone for the fixed carrier density n = 0.75. Fig. 8(a), (b) correspond to the 
cases of a = and a = 1 respectively. 

In case of a = 0, xo, X3 and X5 exhibit the well known incommensurate peaks around the 
nesting vector q = (tt,tt), since these two terms directly reflect the symmetry of the Fermi 
contour. On the other hands X2, and Xi vanish at q = (-7T, tt), as was suggested in §3. The 
reason why Xi also vanish at q = (tt, tt) is that then the integrand of the Xi becomes odd 
function with respect to the energy. 

In case of a = 1, the overall structures and magnitudes of the modified polarizations are 
similar to the case of a = 0. 

It should be noted that in any case Xi and Xa which appear in the numerators of C(q) 
and S(q) in the expressions (|32]) and (|33|) become very small in comparison with the %3 and 
X5 which appear in the denominators. It follows from this observation that we expect the 
terms JCxi and JSx2 to become negligiblly small in comparison with xo- 

Indeed as shown in Fig. 9, JCxi and JSx2 have very small value in the whole region of 
the first Brillouin zone. However we can say with certain that these terms have tendency 
to reduce the magnetic fluctuation, because these terms have negative values in the whole 
region. 

We can see directly from the results shown in Fig. 9 that in case of a = 0, the nesting 
condition strongly surpresses both of JCxi and JSx2- This means the exchange scattering 
processes are strongly surpressed due to the nesting condition. 

On the other hand, in case of a = 1, the exchange scattering is enhanced as the system 
becomes nearer to the half-filling. Furthermore we note that the absolute vale of these terms 
are smaller than that in case of a = 0. We can say this difference comes from the different 
geometry of the Fermi contour. 

Finally we compare Xo and Xo — Xo + JCxi + JSx2- As shown in Fig. 10, we can see 
that we can't distinguish these two quantities. This situation is just what we expect and 
therefore we can certainly say that we can work well only by taking xo into account. This 
means in the present model, the spin fluctuations of the system can fairly well be described 
in terms of the simplified, or effective, RPA, which include only the string type procesesses 
in Fig. 3 in which the single loop can be replaced by the bare one. 

We can expect that this situation can survive in case of the t-J model. In the t-J model 



3 We note that if we descard the vertex correction, it is no longer necessary to include the self-energy 
correction to ensure the Ward-Takahashi identity. 
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the free carrier in the present work is repaced by the free spinon and, in the mean field 
level, the bare hopping parameter t and t' are replaced by t and t' which depend on carrier 
concentration. However fundamental scattering processes are completely the same as in the 
present work. 

In the present paper we have establised the consistent formalism to treat the spin fluctua- 
tions by an itinerant Heisenberg model within the framework of RPA. Here the conservation 
low of the spin and the rotational symmetly of the system was carefully treated. As a result 
of these careful consideration, the simplified version of RPA turns out to work fairly well 
even in the real High-T c problems. 
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Appendix A 

Density of States 

Here we derive an analytic expression for the density of states per spin. 




— I dx I dy5(e + 2t(cosa; + cosy — a cos x cosy)) 

IT 2 JO JO 



2vr 2 J-i ^(1 - £ 2 ){(1 + e/2) + (1 - a)£}{(l - 7/2) - (1 + a)£} 



1 K[ 4-(a-e/2t)2 
7T 2 v/2+^ [ \ 2(2 + as) 1 



where K(x) denotes the elliptic integral of the first kind. Here we made use of a formula [|12|; 




2 



F[sin 



(a — c) (b — u) {b — c) (a — d) 




\ (b — c) (a — u) ' \ (a — c) (b — d) 
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and the relation F(7r/2,g) = K(x) where F(k, x) denotes the elliptic integral of the second 
kind. Here a>b>u>c>d. 



Appendix B 

Bethe-Salpeter Equation for particle-hole Channel 

Here we present the derivation of the solution of the Bethe-Salpeter equation for particle-hole 
channel fl3"U|). The equation to be solved is 

T(k + q,k) = l + ~Yl Jv-vG{k' + q)T(k' + q, k')Q{k'). (52) 

k' 

By noting that 

Jk-k' = 2 J(cos k x a cos k' x a + sin k' x a sin k' x a + cos k y a cos k' y a + sin k' y a sin k' y a), (53) 
we can rewrite (0) in the form, 

T(k + q, k) = 1 + JC(q) (cos k x a + cos k y a) + JS(q) (sin k x a + sin k y a). (54) 

Here 

1 1 

-(cos k' x a + cos k' y a) , S(q) = - 

where 



C(q) = -(cos k' x a + cos k' y a) , S(q) = -(cos k' x a + cos k' y a), (55) 



^^^^(fc' + g,^^'), (56) 

JV fc' 



and C(g) and 5(g) are to be determined self-consistently. By inserting (|55f) back into (|52|). 
we obtain the following equations, 

C(q) = -\xi{q)-\jC{q)xM-\jS(q)xM), (57) 

S(q) = -\x2{q)-\jC{q)xM-\jS{q)X^q). (58) 

where Xii ' ' ' > Xs have been defined in (|3"4T). Finally by solving these simultaneous equations, 
we obtain (R2) and (B3). 



14 



References 

[1] H.Won, and K.Maki, Phys. Rev.B.49(1994), 15305. 

[2] P.W.Anderson, Science 235 (1987), 1197; 

F. C.Zhang and T.M. Rice, Phys. Rev.B.37 (1988), 3754 

[3] G.Baskaran, Z.Zou, and P.W.Anderson, Solid State Communications, 63, 973(1987). 

[4] Tanamoto.T,Kohno.H,Fukuyama.H,J.Phys.Soc.Jpn.60(1991),3072; 
Tanamoto.T,Kohno.H,Fukuyama.H,J.Phys.Soc.Jpn.62(1993),717; 
Tanamoto.T,Kohno.H,Fukuyama.H,J.Phys.Soc.Jpn.63(1994),2739. 

[5] P.W.Anderson, Phys. Rev.ll0(1958),827 

[6] N.F.Berk, and J.R.Schrieffer, Phys. Rev. Lett.l7(1966),433; 
S.Doniach, and S.Engelsberg, Phys. Rev. Lett. 17(1966), 750; 
P.W.Anderson, and W.F.Brinkmann Phys. Rev. Lett. 30(1973), 1108; 
S.Nakajima, Prog. Theor. Phys.50(1973),1101. 

[7] P.A.Wolff, Phys. Rev.l20(1960), 814; and as a review e.g. 

G. D.Mahan, Many- Particle Physics, Sec. 5. 8, (Plenum Press, 1990) 

[8] A.A.Abrikosov, J.C.Campzano,K.Gofron,PhysicaC214(1993),73 

[9] J.Yu and A.J.Freeman, J.Chem.Solids 52(1991), 1351 
[10] B.S.Shastry, J.Phys.F,Metal Phys.9(1979),1367 
[11] G.Baym and L.P.Kadanoff, Phys. Rev.l24(1961), 287; 

[12] I.S.Gradshteyn and I. M. Ryzhik, Ta ble of Integrals, Series, and Products,Ac&demic Press, 
Inc.1980. 

Figure Captions 

Fig. 1: From left to right, the energy contour, corresponding DOS, and the structure of the 
saddle point near the Y-point for 

(a) a = 0, 

(b) a = 0.5, 

(c) a = 1. The thick line in each figure represents the Fermi contour at the half filling and 
the broken line corresponds to the energy level at which the DOS becomes divergent. 
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Fig. 2: The fundamental processes induced by J q . 

(a) The scattering process with spin flip (type-a), 

(b) the scattering processes between parallel spins (type-b) and 

(c) the scattering processes between anti-parallel spins (type-c). Here the straight line and 
the wavy line represent respectively the green's function of an itinerant carrier and the 
antiferromagnetic interaction. 

Fig. 3: The string type summation for the transverse susceptibility. 

Fig. 4: Bethe-Salpeter equation for the triangle vertex. 

Fig. 5: The process corresponding to XoOn)- 

Fig. 6: The process corresponding to ( fiTf ) in the text. 

Fig. 7: Temperature dependence of X(T) in the self energy for 

(a) a = 0; n = 0.85 and 

(b) a = 1; n = 0.85. We set J = t/4. 

Fig. 8: Scan of Xo, '''5X5 along the r — > X — > M — > T line in the first Brillouin zone for 

(a) a = 0; n = 0.75 and 

(b) a = 1; n = 0.75. We set T = O.OAt ~ 60 [K] 

Fig. 9: JC(q)xi(q) and JS(q)x2{.<l) in the first Brillouin zone for 

(a) a = 0; n = 0.95,0.75 and 

(b) a = 1; n = 0.95,0.75. 

Fig. 10: Xo(?) an d i n the first Brillouin zone for 

(a) a = 0; n = 0.95, 0.75 and 

(b) a= 1; n = 0.95,0.75. 
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